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While general relativity is an extremely robust theory to describe the gravitational interaction in 
our Universe, it is expected to fail close to singularities like the cosmological ones. On the other 
hand, it is well known that some dark energy models might induce future singularities; this can 
be the case for example within the setup of the Holographic Ricci Dark Energy model (HRDE). 

On this work, we perform a cosmological quantisation of the HRDE model and obtain under which 
conditions a cosmic doomsday can be avoided within the quantum realm. We show as well that this 
quantum model not only avoid future singularities but also the past Big Bang. 
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I. INTRODUCTION 

The late-time observed acceleration of the Universe has promoted several scenarios that try to describe this recent 
speed up of the Universe. In fact, it is well known that a Universe filled only with matter and radiation (for open, 
flat or closed spatial geometries) cannot expand with positive acceleration. Therefore, it is necessary to find other 
mechanisms/matters that could explain this feature. The simplest phenomenological approach consists in invoking 
another kind of energy density responsible for the current acceleration of the Universe which is usually dubbed dark 
energy (DE) being the cosmological constant the simplest option [1|. 

The common way to describe DE is via its equation of state (EoS) parameter, which is usually denoted by ude- 
This cosmological parameter is the ratio between the pressure and the energy density of DE. It can be constant or 
time-dependent. By definition, DE EoS must fulfil wde < —1/3 at late-time to be able to describe the current speed 
up of the Universe. In fact, observations suggest wde ~ — 1 at present (see for example i)- 

Within the above framework, the ACDM model is the best fit to the observational data 0. This model assumes 
besides ordinary matter, the presence of non-baryonic matter corresponding to dark matter (DM) and a cosmological 
constant A. In fact, the EoS of A is just wa = —1 and remains unchanged in time Q. In this Universe, the expansion 
is accelerated at late-time reaching an infinite scale factor in an infinite cosmic time where the Universe geometry 
is described by a de Sitter space-time. Planck latest results when combined with other cosmological measurements 
provides wde ~ —1.006 ± 0.045 for a DE constant equation of state @. Even though this value is very close to minus 
one, a slight deviation from it is crucial for the asymptotic future behaviour of the Universe. In fact, if wa is slightly 
larger than -1, the Universe will continue expanding eternally, but if this value is slightly smaller than -1, and no 
matter how tiny is such a deviation, the Universe could end up in a doomsday. For example, it could happen that 
the energy density becomes infinite at a finite cosmic time I4H12||. 

Although the ACDM model gives the best observational fit to explain our current Universe, there are other ways 
that have gained great attention [I]. One of them is the holographic dark energy scenario mill which is inspired 
on the holographic principle rooted in quantum gravity. We next explain briefly the ideas behind the HRDE. 

As it is well known, the entropy of a given closed system with finite volume L 3 has an upper bound which is not 
proportional to its volume, but to its surface area, L 2 0®. On the other hand, for an effective quantum field 
theory with a given ultra-violet (UV) cutoff, Muv, the entropy of the same system is scaled as L 3 M^ V . Therefore, 
there is always a scale or a length where the quantum field theory with UV cutoff is expected to fail. This will happen 
for large volumes or lengths. To overcome this problem a link between UV and infrared (IR) cutoffs was proposed in 
0 : 


L 3 m£ v < LM 2 p . 


( 1 . 1 ) 


This proposal ensures the validity of the quantum field theory within this regime. When the inequality is saturated, 
we can define an energy density which is inversely proportional to the square of the characteristic length of the system. 
Applying these ideas to the universe give rise to the holographic dark energy scenario ®. 

Now, the next task would be to find a suitable holographic energy density able to speed up the current Universe. 
This issue has been tackled in different works. For example, taking as a characteristic length the inverse of the Hubble 
parameter, In = i? _1 , it was found that the effective EoS for DE is equal to zero and therefore is not a suitable 
proposal to describe the current Universe because it would imply an eternally decelerating Universe [H}. Later 
on, the particle horizon was suggested as a characteristic length for the universe within the holographic approach 
Zph = a fg dt/a ®. In this case, it turns out that the effective EoS is larger than -1/3, therefore this choice is equally 
unsuitable to describe our present Universe ® ®- Contrary to the previous proposals, the future event horizon 
Z E h = a, J* f dt/a is phenomelogically viable as it fits the current observations (l3Hl9|. however, it has a drawback with 
causality; the future event horizon should not affect the current or past physical evolution of our Universe [0. There 
is a further possibility to define an holographic dark energy model which consists in taking as the square of the length 
characterising the Universe, the inverse of the Ricci scalar curvature [21]. This model was named the Holographic 
Ricci Dark Energy scenario (HRDE) for further generalization of the holographic DE model see ]22 

While the HRDE is suitable to describe the current acceleration of the universe as shown in [2ll . 
induce a Big Rip singularity [2lj| : i.e., the scale factor, the Hubble parameter and its first cosmic time derivate blow up 
in a finite future cosmological time • This model has been observationally constrained in [33]. For more updated 
observational constraints on the HRDE we refer to Ref. where even interaction between DM and the HRDE is 
considered, and again, a Big Rip is favoured observationally. 

When the Universe approach the Big Rip regime, we expect quantum effects to be important, therefore, it is 
necessary to do a quantum analysis. The field of cosmological singularities has been extensively studied in quantum 
cosmology [34], ®] (see [36] for a review on this topic). Quantum Cosmology consists in applying the quantum 
theory to the Universe as a whole [13] ■ A consistent theory of quantum gravity should in principle avoid the classical 



it might 
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singularities prevalent in the classical theory of general relativity [38l440|. Indeed, it has been shown that in some 
specific models, most DE related singularities can be avoided in the analogous quantum version [38l - |4ll ]. 

In the present work, we address the quantisation of the HRDE model. We use the Wheeler-DeWitt (WDW) 
formalism for a homogenous, isotropic and spatially flat universe. The solutions to the WDW equation must obey 
the DeWitt boundary conditions [34] . i 35| . which implies that the wave function of the Universe has to vanish close 
to the singularities, ensuring that the classical singularity is avoided through the quantisation procedure. We extend 
our analysis to the primordial Universe where a Big Bang is expected to take place from a classical point of view and 
where again a quantum analysis is required. 

The paper is outlined as follows. In the next section, we make a brief review of the HRDE model. In section III, 
we present the WDW equation for a standard fluid and present the solutions for the HRDE model. Finally, in section 
IV we discuss the overall conclusions. In the Appendix, we include a brief explanation of the WKB method that was 
used to solve the Wheeler-DeWitt equation for some periods of the expansion of the Universe. 


II. THE HRDE: A SHORT REVIEW 


We start reviewing the HRDE model. The Universe on its largest scale can be described by a FLRW Universe 
whose metric reads 


ds 2 = —dt 2 + a 2 (t) 


dr 2 


+ r 2 d0 2 + r 2 sin 2 9dip ) , 


1 — kr 2 


( 2 . 1 ) 


where k = —1, 0,1 for open, flat and closed spatial geometries, respectively. Aside from the HRDE, our Universe will 
be filled with cold dark matter (CDM), baryonic matter and radiation. Therefore, the Friedmann equation reads 


H = 


8t tG 


E fC 

Pi - 7^’ 


( 2 . 2 ) 


where H = a/a is the Hubble parameter and dot denotes a derivate with respect to the cosmic time. The summation 
runs over radiation, matter, and HRDE. The energy densities of radiation and matter can be described as follows 


_3i?o 2 o ( a 

Pr ~ SttG r ° Uo 


-4 


(2.3) 


3 H, 


Pm — 


_0_Q 

87rG \ ao 


-3 


(2.4) 


where Q r o and fl m o are the dimensionless energy density parameters at present time for radiation and the matter 
components (cold dark matter and baryonic matter), respectively. In turn, Ha is the value of the Hubble parameter 
at present. For the HRDE, the expression for the energy density is the following [21 ] 


p R = 6/3 [H + 2H 2 + 


(2.5) 


where f3 is a positive constant. A dot stands for derivates with respect to the cosmic time t. Defining a dimensionless 
quantity J3 = 167rG/3, and solving the Friedmann equation, the expression for the Ricci dark energy density is found 
to be [2l[ 


PR 


3Hl 

8t rG 



( 2 . 6 ) 


where fi p o is an integration constant which will quantify the effective amount of DE in the HRDE model [2l| . Notice 
that the presence of radiation and spatial curvature in the Universe do not modify the previous result. The Ricci dark 
energy has one part that behaves as matter and another part which depends on the value of /3 and plays the role of 
DE. The asymptotic future behaviour of the Universe depends on the values acquired by /3, more precisely: 


1. If 1 < /? , the cosmic acceleration is negative. We disregard this case as it cannot describe the present Universe. 
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2. If 1/2 < /3 < 1, the Universe enters in an accelerating state when the HRDE dominates. The Universe is 
asymptotically flat in the future. 

3. If /? = 1/2 , the model is equivalent to the existence of a cosmological constant plus the matter contributions. 
We will disregard this case as it reduces to ACDM model. 

4. If 0 < /3 < 1/2, the Universe not only enters in an accelerated state, but also super accelerates (H > 0) in the 
future hitting a Big Rip; i.e. the Universe hit a singularity at a finite cosmic time. 

The observational constraints of the HRDE model favour the last case, i.e 0 < /? < 1/2 [32, ; 33|, so the Universe 
would reach a future singularity in a finite time. In this framework, classical Einstein theory is no longer valid and it 
is necessary to make a quantum treatment. 


III. QUANTISATION OF THE HRDE 


The Wheeler-DeWitt equation is the analogous to the Schrodinger equation taking the Universe as a whole; i.e. as 
the system to be analysed [H HU- In order to deduce the WDW equation, the first step is to obtain the classical 
Hamiltonian from the usual Hilbert-Einstein action (in our case we restrict to a FLRW metric). Then the quantum 
operators are introduced according to the canonical quantisation procedure leading to the WDW equation [HI, HU. 
The gravitational action depends upon the chosen metric, in particular for a FLRW space-time, the Hilbert-Einstein 
action is given by 


She = 


167tG 


d 4 a 


■y/=g(R- 2A)-2 


d 3 x\/—hK 


(3.1) 


The second term on the right hand side of the previous equation contains the extrinsic curvature K and is a compulsory 
boundary term to have a well-defined variational problem. The extrinsic curvature and its trace read 


K ab = 


1 dh ab 
2 N dt ’ 


(3.2) 


K = K ab h ab , 

where N (t) is the lapse function and h ab is the induced spatial metric 

dr 2 


hab = a (t) 


1 — kr 2 


+ r d6 + r shr 9 dip 


(3.3) 


(3.4) 


The Lambda function introduced in the action corresponds to A = 8tt Gp\ i.e. the matter Lagrangian of the Universe 
is described by A (a). From the Hilbert-Einstein action, we obtain the following Lagrangian density (after performing 
an integral) 


L = N 


3n ( ad 2 , . , , a 3 

_(__ + fco _A(a) T 


Defining 


dL 

Pa — TT = 
oa 


37t f aa\ 


2G \ N J ’ 


we get the Hamiltonian corresponding to the action m 

u = n 


g_pI 

37t a 


-^ + ^A(«)«3 


(3.5) 


(3.6) 


(3.7) 


where p a is the canonical momentum corresponding to the scale factor. Notice that the canonical momentum associated 
with the lapse function does not appear because it constitutes a primary constraint (35|- Therefore, the scale factor 
is the only variable in our cosmological problem. The curvature term is reabsorbed in the A(a) function; i.e.: 


A (a) = 3 H 2 


n 


r 0 


ao 


-4 


+ 


Q 0 (— 

2-/3 ) m °Uo 


-3 


+ fife' 


fcO 


ao 


a 


pO 


ag 


~ 2 ( 2 -})' 


(3.8) 
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Here, Qko is the dimensionless energy density parameter for curvature at present. 
The term ]r a /a generates the operator 



a 


-hf 


a *d a 


a id a 


(3.9) 


in the quantum framework where we have choosen a factor ordering corresponding to the covariant generalization of 
the Laplace-Beltrami operator ||34j (for alternative choices see for example (42l - l44l |). It is useful to apply the following 
change of variable to remove the first order derivate from the quantum Hamiltonian operator 


x = 



(3.10) 


Therefore, the quantum Hamiltonian operator can be written as 


H = N 


f 3 Gh 2 2 SitHqCLq 


i 4 ™o' 


AG 


Q r QX 3 H - 


2-/3 


T ^kQX 3 T D^qX 3 ^ 1 ^ 


(3.11) 


As the variation of the Hamiltonian with respect to the lapse function N produces the Hamiltonian constraint, the 
WDW equation reads 


HA/(x) = 0. (3.12) 

We will take the case of a spatially flat universe (flko = 0) for simplicity and in accordance with the current observations 
[2]. Therefore the WDW equation reduces to 


dl + 7 


fl. 


qX 3 + 


2-/3 


H m0 + n p0 x 


T (x) = 0 , 


(3.13) 


where we have introduced the following dimensionless constant: 

7 = (nHoal/Gh) 2 . (3.14) 

Due to the complexity of the equation (13.131) . we will divide the evolution of the Universe thereof in three regimes 
corresponding to the domination eras of radiation, matter and DE, respectively. The first one, has an exact solution, 
for the others two cases we will make a WKB approximation (up to first order). In the radiation dominated era the 
WDW is 


{ d 2 x + 7U r0 x 3 1 'k (;r) = 0, 

whose exact solution reads [45| 


Ti(x) =(D r 0 7 ) 8 \[x 


C1J1 -\AW72; 3 +C2Y3 


rOjX- 


(3.15) 


(3.16) 


where C\ and C 2 are constants. The functions J 3/4 and I 3/4 correspond to the first and second kind Bessel functions 
of order 3/4, respectively. We choose C 2 = 0 to ensure that the wave function vanishes when a —> 0, according with 
the DeWitt boundary condition [33;, 34|. For the matter dominated era the WDW is 

{dl+ig 2 {.x)}A> (x) =0, (3.17) 


where the function g 2 {x) is defined as 


92 (x) 


Q r QX 3 + 



(3.18) 


The first order WKB solution gives 


A>2 (x) « [-732 ( a ;)] 4 


ai e ih2{x) + a 2 e~ ih2{x) 


(3.19) 
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where <37 and a 2 are constants and the function h 2 {x) is given by 


h 2 (x) = 4/7 


2-/3 


Cl. 


mO 


X3 


2 - B\ Cl 


2 Cl 


L r 0 


^mO 


Finally, during the DE dominated era the WDW is; 

+753 (x)} 4- (x) = 0, 

where the function gs(x) is defined as 

2 


53 (x) = 


2-/3 


^mO + ^pO X 


The first order WKB approximation gives the solution 

*3 (X) « [-753 (®)] _i [<5ie iM:c) + S 2 e~ lh ^ x) 
where S\ and S 2 are constants. The function hz{x) reads 


h * (x) = ^x | 2 v ^M+ a /(^) 


^mO 1 


1 1 1 f &- 2 \ a 

77 ’ ? i- T — 5 


2 s s V 2 7 fi. 


^pO 


£ m 0 


(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 


where s = —2/3(1 — 2//3). Notice that (i) the hypergeometric function defined in (13.2411 is not well defined as a series 
but it is well defined as an integral [45| and (ii) the parameter s is positive, in fact it is larger than 2 for 0 < /3 < 1/2. 
The function h%(x) is always real, on the other hand, the function ( 73 (x) is positive and an increasing function of 
x\ in fact, it blows up for large values of x. Therefore, the wave-function ()3.23ll which corresponds to an Universe 
dominated by DE vanishes at large scales and the DeWitt condition is fullfilled authomatically for 4 / 3 . 

Now, it is necessary to connect continuously the solutions among them. Taking the arbitrary constant C7 = 1 
(for convenience, this election will not modify the fundamental behaviour of the wave function), the conditions for a 
smooth wave-function give the values of 07 , a 2 , Si, S 2 . These conditions are just the continuity conditions of the wave 
function and its first derivate on a first connecting point ( 27 ) and on a second connecting point (X 2 ), which can be 
read as 


Ti ( 27 ) = 4 / 2 (xi), 4>i/ ( 27 ) = 4V (xi), 4 > 2 (x 2 ) = 4 > 3 ( 27 ), 4V (x 2 ) = 4v ( 27 ), (3.25) 


where prime stands for a derivative with respect to x. 

We choose the first connecting point ( 27 ) in which the matter component is subdominant with respect to the 
radiation component (for example p m ~ 10 ~ 4 p r ), the second connecting point ( 27 ) is just where the radiation and 
the phantom contribution in equation (13.131) are subdominant with respect to dark matter (for example when they 
are equal). The connecting points then read 


Xl 


\( 2 ) 

^mO 

K2-/3J 

^rO 


3 

2 


(3.26) 


x 2 = 


M 


Using the Cramer method to solve the algebraic system, the constants 07 , a 2 , <5i, S 2 can be written as follows 


on = — 2 [-752 (xi)p 


V 1 ai 2 
52 a 22 


(3.27) 


(3.28) 


ot 2 = 


2 [ — 752 (xi)] 2 


Sr =- 


2 [-753 (x 2 )]< 


an y 1 
a 2 i 52 

Zl bi2 

z 2 b 22 


(3.29) 

(3.30) 





































where we define 


^2 = ~2 [-733 (* 2 )] : 


6n zi 
&21 -22 


an = [-732 (*i)] _i e ih2 ^\ 
ai2 = [-732 (si)H e~ lh ^ xi \ 


7 

(3.31) 

(3.32) 

(3.33) 


« 2 i = | [- 32 (a:i)] 4 (|firoaa 3 ) -*[ 732 ( 37 )] 2 \ e 


l J h 2 (^ 1 ) 


(3.34) 


a 2 2 = < [-32 ( 27 )] 1 3 ) +i [ 752 (^ 1 )] 2 } e 


“I l p -ih 2 {xi) 


(3.35) 


31 = (^ 07 ) 8 


(3.36) 


32 = (^r07) 8 


2 ^ 


Ja ( 2 \/^ro7^i ) - 2 (^ 7 - 07 ) 2 Xi 3 Jt, ( - V ^07^1 


+ 1 


(3.37) 


6 n = [-733 (* 2 )]" 4 e ift3 ^ 2 ), 
6 12 = [-'yg 3 (x 2 )]-^e- lh3 ^\ 


(3.38) 

(3.39) 


6 2 1 = 


J 7" 


[- 33 (^ 2 )] 4 (sn p0 x s 2 X ) - i [733 (x 2 )] 2 l e ^ 3 (a = 2) 5 


(3.40) 


&22 = V [-53 (^r 5 (sO^r 1 ) +ihg 3 (x 2 )]-He- ih3 ^\ 


(3.41) 


Zi = [-732 (x 2 )]' 


ai e ih ^) + a 2 e- ih2 ( X2) 


(3.42) 


^2 = ai < [-32 (x 2 )] 4 Qfiroa: 2 3 ) - * [732 (x 2 )] 2 |> e lh2< - X2) 


-a 2 | [-32 (x 2 )[ 4 ^^0*2 3 ) + * [732 (z 2 )] 2 } e 


— 3 l p -ih 2 {x 2 ) 


(3.43) 

(3.44) 


With the aim to show a numerical result, we use the constraints /3 = 0.3823, fl r o = 8-10 , fl m o = 0.2927, fl p 0 = 0.6380 

in accordance with the best fit of the current Universe within the HRDE model [33]. Therefore, the connecting points 
are 


x 1 = 3.287 • 10" 12 , x 2 = 0.07608, 


(3.45) 


and the constants can be written as 


07 — Uq T UqZ, q?2 — Uq WqZ, $r — u?j H - **,5*, S2 — ~t - u^i. 


(3.46) 


















We next show an example of these constants. We choose 7 = 10 20 , sufficiently large to ensure the validity of the 
WKB approximation for the connecting first point; i.e. well inside the radiation dominated epoch where matter is 
subdominant. Notice that if the WKB approximation is fulfilled for the first connecting point it is well defined for 
the second connecting point. Then we obtain q(x\) = 0.035 (see Eq. (|A4[) h Finally, using the above introduced 
constraints and for the selected 7 we find: 

u a = 80274.39, v a = -288039.58, u s = -298682.24, v s = 14131.34. (3.47) 

In principle, 7 is much larger, in fact, taking ao = 1 gives the value 7 = 8.072 • 10 8 '. If in addition, we take into 
account the whole number of e-folds since the radiation dominated epoch, the parameter 7 would be much larger. We 
choose the previous value of 7 ; i.e. 7 = 10 20 , to be able to get numerically the values of the constants defined in Eq. 

(G 3 SD■ 


IV. CONCLUSIONS AND OUTLOOK 

The HRDE is a suitable proposal to describe the late Universe. The best fit of the model provides a value for the 
proportionality constant /3 inside the interval 0 < /3 < 1/2 [33|. This means that the Universe is not only accelerating 
but will also face a future Big Rip singularity. This classical singularity is analysed within a quantum treatment 
where the quantisation is realised in the framework of the WDW equation for a flat FLRW universe and imposing 
the DeWitt boundary condition. We have shown that the Big Rip singularity as a consequence of the phantom like 
behaviour of the HRDE could be avoidable and would be harmless within the quantum approach used in this work. 
This might not be the case within a classical approach (46) . The DW condition can be regarded as a guidance in 
the nowadays incomplete theory of quantum cosmology. In fact, the disappearance of the probability distribution at 
singularities should arise in the theory in a natural way as a dynamical consequence of some other requirements, such 
as the normalizability of the wave function, and should not be postulated. Given that we lack of a complete and 
consistent quantum gravity theory, we will stick to the DW condition as our guidance for singularity avoidance. 

Despite the quantum analysis presented in this work about the avoidance of the Big Bang and Big Rip singularities, 
this fact cannot be interpreted as an exact and thorough evidence of singularities avoidance in quantum cosmology, 
but rather an indication that a consistent and complete quantum theory of gravity should be free of these singularities. 

To solve the Wheeler-DeWitt equation we carry two types of approximations: (i) we divide the evolution of the 
universe in three different epoch corresponding to radiation, matter and DE dominance as explained in the previous 
section, (ii) for the last two periods we use a WKB approximation where it is enough to go to first order in the 
WKB approximation to ensure the DeWitt boundary condition for large scale factors where the wave function is 
asymptotically decreasing and vanishing. 

During the radiation dominated epoch, we obtain the exact wave function, fulfilling the DeWitt condition, which 
can be matched with the WKB approximation for the second period (matter+radiation). The larger is the value of 7 , 
defined in (13.141) . the sooner, i.e. for smaller scale factors, we can carry the matching between the two wave functions. 
In fact, the asymptotic behaviour of the wave function during the radiation dominated epoch matches smoothly and 
naturally with the WKB approximated solution corresponding to the matter and radiation epochs. 

The quantisation is necessary to describe the Universe close to singularities. However, outside of these singularities 
the Universe can be described classically. In this regime, the square of the modulus of the wave function can be 
interpreted as the classical probability density distribution. We can see this clearly when the value of 7 is very large 
and therefore, the wave function carries out a significant number of oscillations within a short interval of the chosen 
variable. In fact, these oscillations are modulated by the classical probability density distribution, which is defined as 
the time average of the scale factor. 

In an analogy with the classical picture, a slight deviation of the resulted wave function from the classical probability 
density distribution is expected. This is due to the performed first type of approximation, which disregards the 
contribution of DE for small scales and the contribution of radiation for large scales. Therefore, this deviation is more 
significant in the matter dominated epoch. 

Next, we would like to stress that our model has only one degree of freedom described through the scale factor 
which in fact can play the role of the classical time. On the other hand, gravity is a reparametrisation invariant theory 
with first class constraints. Therefore, there is always a gauge fixing condition which reduces the number of physical 
variables and in our case we would be left without any degree of freedom after the gauge fixing. However, our current 
work can be regraded as a first approach in quantising the HRDE model, in fact a toy model, and we will present 
in a different work a more elaborated scenario for the quantisation of the HRDE with two physical variables and a 
genuine degree of freedom }47j ] . In fact, to get at least one degree of freedom with physical meaning, we can map the 
matter content given by a perfect fluid to one or more scalar fields that could mimic the different components of the 
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universe. This method was carried out for example in [48| for a minimally coupled scalar field or a tachyon scalar 
field. 

In addition, we have been mainly focussing on the wave function of the universe while in fact the important thing 
is the probabilistic interpretation of it. This requires the definition of a Hilbert space with a proper scalar product 
and measure. In the interesting review [4£j, this non trivial problem is discussed and several potential solutions are 
presented while in mm those procedures are applied to different cosmological models. We hope to implement those 
methods for the HRDE in a future work. 

The aforementioned interpretation has some conceptual drawbacks: first, this picture corresponds with the descrip¬ 
tion of an “external observer”, who lives in a reality outside from the quantum system under study. This is certainly 
not the case for a cosmological system where the observer is part of the system |49[. In fact, we would need to apply 
a quantum theory of closed systems or the many-worlds interpretation of quantum mechanics to our universe but this 
is beyond the scope of our current work (see [49( for more details on this subject). 

We have restricted our analysis to a homogeneous and isotropic configuration. However, it is well known that the 
inhomogeneities and anisotropies can be quite important close to singularities, for example through the creation of 
particle and gravitons (see for example [Kll, [52])- We leave this interesting issue for a future work. 
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Appendix A: Justification of the WKB approximation 

For the following differential equation 

{dl+19 0 )} ^ 0 ) = 

the first order WKB approximation gives the solution [53] 

T (x) » (<? (x))“* [Aie^ + A 2 e~ h( - x) 1 , 
where A\ and A 2 are constants of integration and h{x) reads 

M*) = /VAFM< 1 *. 

The first order WKB approximation is valid in the region that complies with the following inequality 


q(x) = - 
7 


5[gf (x)] 2 3 4 — 4 [gtt (x)] [g (x)] 


16 [g (x )]' 1 


< 1 . 


(Al) 


(A2) 


(A3) 


(A4) 


For the chosen interval (13.45)) . this inequality is fulfilled for large values of the dimensionless quantity 7 . On the 
previous expression a prime stands for a derivative with respect to x. 
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